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Abstract 

Millimeter wave (mm-wave) and massive MIMO have been proposed for next generation wireless 
systems. However, there are many open problems for the implementation of those technologies. In 
particular, beamforming is necessary in mm-wave systems in order to counter high propagation losses. 
However, conventional beamsteering is not always appropriate in rich scattering multipath channels with 
frequency selective fading, such as those found in indoor environments. In this context, time-reversal 
(TR) is considered a promising beamforming technique for such mm-wave massive MIMO systems. In 
this paper, we analyze a baseband TR beamforming system for mm-wave multi-user massive MIMO. 
We verify that, as the number of antennas increases, TR yields good equalization and interference 
mitigation properties, but inter-user interference (IUI) remains a main impairment. Thus, we propose a 
novel technique called interference-nulling TR (INTR) to minimize IUI. We evaluate numerically the 
performance of INTR and compare it with conventional TR and equalized TR beamforming. We use 
a 60 GHz MIMO channel model with spatial correlation based on the IEEE 802.llad SISO NLoS 
model. We demonstrate that INTR outperforms conventional TR with respect to average BER per 
user and achievable sum rate under diverse conditions, providing both diversity and multiplexing gains 
simultaneously. 
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I. Introduction 

Massive MIMO systems have been recently recognized as one of the technologies that can 
bring unprecedented performance gains for next generation wireless communications [jT|. Among 
its potential benefits, noise, fading and inter-user interference (IUI) effects have been shown 
to progressively reduce as the number of antennas in the system increases [J2]|. Thus, a large 
number of antennas simplifies the multiple access layer and increases the system’s capacity [0. 
However, many challenges remain for the implementation of massive MIMO systems such as: 
channel estimation and reciprocity issues, large pilot overheads, hardware cost, size and power 
limitations, network architecture adaptations, antennas and propagation aspects 0. 

Recently, millimeter wave (mm-wave) and massive MIMO have been proposed in tandem for 
next generation systems 0,0. This can be easily justified because a large number of antennas 
operating at mm-wave frequencies (e.g. 28, 38, 60 and 73 GHz) can be used in compact devices 
due to the small wavelength (4 to 10 mm approx.) and (hence) small antenna sizes. In addition, 
it has been shown that mm-wave networks are suitable for dense small cells (especially in indoor 
environments), as inter-cell interference is naturally mitigated due to high propagation losses at 
those frequencies 01, 0. Another benefit of using mm-wave is the huge bandwidth availability, 
with some standards planning to operate with more than 2 GHz bandwidth, e.g. ||9l . 

Nevertheless, a number of problems arise for mm-wave massive MIMO systems. In particular, 
their performance is highly dependent on the antenna array configuration and the propagation 
environment. Hence, factors such as the coupling between antennas and channel spatial corre¬ 
lation play a significant role on the actual capacity and diversity gain that this kind of systems 
can achieve 0. 

Beamforming is needed in mm-wave systems because of large propagation losses. Tradi¬ 
tionally, antenna arrays use beamsteering techniques in order to increase the received power 
in specific directions iflOll and, consequently, achieve diversity gain. This is usually performed 
with either analog (RF) or digital (baseband) phase shifters in each antenna. Recent approaches 
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to beam forming^] in mm-wave massive MIMO use hybrid analog beamsteering combined with 
digital precoding techniques assuming narrowband fading channels 0, lfl4l - llT7l . This hybrid 
analog/digital solutions are necessary given that fully digital solutions requires one digital to 
analog converter per antenna, which is extremely constly in terms of power. However, con¬ 
ventional beamsteering is not always appropriate in multipath channels with frequency selective 
fading, such as those found in indoor environments. In those cases, more sophisticated techniques 
are required to take full advantage of the number of elements in the array and also multipath 
propagation. 

The specific propagation aspects of mm-wave systems have been recently studied. Statistical 
models for mm-wave channels have been developed in Ifl8ll . Ifl9l . where extremely narrow 
antenna radiation patterns are considered using massive MIMO. These models provide charac¬ 
terization of scattering clusters in the angular and delay domains, power-delay profiles (PDPs), 
and propagation losses in outdoor scenarios. Similar models can be found for indoor scenarios 
lt20l . Another popular SISO channel model for indoor mm-wave systems is the IEEE 802.1 lad 
EO, which considers extremely narrow radiation patterns and analog beamsteering. However, 
there are only few studies on the spatial correlation in mm-wave MIMO channels. An interesting 
work is |[22l . where it is demonstrated that correlation at 60 GHz can be very high due to the 
small number of multipath components (MPCs). It has also been recognized that the specific 
structure of spatial correlation is highly dependent on the scattering environment. Given this 
conclusions, it is not clear yet whether diversity and/or multiplexing schemes should be used in 
order to maximize the system’s gain lf23l . 

In this paper we propose a novel time-reversal (TR) based ||24]| solution to the multi-user 
beamforming problem for indoor scenarios in mm-wave massive MIMO, which provides both 
diversity and multiplexing gains. TR is a transmission technique that enables spatial focusing of 
the signal at the receiver by using the time-reversed channel impulse response (CIR) as a linear 
filter applied to the transmitted signal llT2ll . lH3ll . TR is considered a promising technique for 
future massive MIMO systems @. By using TR, all multipath components are added in phase at 


'The term beamforming is traditionally used to denote phased array techniques for beam steering, i.e. operating in the 2D 
manifold spanned by the azimuth and elevation angles. In this paper, we shall use the term beamforming in a broader sense 
to denote signal processing techniques that allow spatial focusing of RF power in co-range as well (3D) or even in time (4D 
space-time beamforming) rm-fm 
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the receiver at a specific instant providing i) an increase in the signal power in the surroundings 
of the receiver (commonly referred to as spatial focusing), and ii) a partial equalization effect 
(commonly known as time focusing) that reduces inter-symbol interference (ISI) caused by the 
channel’s frequency selectivity [|25l . If26l . This features enable low computational complexity 
receivers, which is a key advantage of TR with respect to multicarrier (OFDM-like) systems ll27l . 
Moreover, multipath components add incoherently at regions in space away from the receiver, 
mitigating interference to other users f28j|, lf29ll . 

A number of works have addressed different aspects of TR beamforming, with particular focus 
on single user systems If26l . If28l - lf30ll . In these references, the spatial and temporal focusing 
properties of TR have been considered, and both theoretical and empirical characterizations of 
bit error rate (BER) have been made under specific scenarios and channel models. A common 
finding in the literature is that ISI is the main limiting factor of TR. This is because ISI imposes 
a lower bound in the achievable BER at high signal to noise ratios (SNR) in single user systems 

m. 

The challenge of mitigating ISI in TR has also received increasing attention. Different equaliz¬ 
ing solutions have been proposed in ll29ll . OTft . Il32l for single-user systems. An important result 
of |[29l is that the ratio between the desired signal power and the ISI power in TR increases 
linearly with the number of antennas. Thus, BER performance can potentially have a significant 
improvement when TR is applied in massive MIMO, without additional equalization. 

For multiuser systems, TR for multiple access in the downlink was proposed in ll33l and Il34l . 
where IUI is recognized as the main limiting factor of BER performance. Also, lt24l proposes 
several multi-user TR techniques. A multiple access TR technique that uses rate-backoff is 
proposed in lf35l . where an approximation for the signal-to-interference-plus-noise-ratio (SINR) 
is given, showing that it increases with the number of antennas. 

However, previous works have not addressed the following aspects: 

. Proposed beamsteering techniques in mm-wave massive MIMO are narrowband (for flat 
fading channels), and do not take take full advantage of multipath propagation to increase 
diversity gain. Thus, these techniques may not be appropriate for frequency selective chan¬ 
nels found in indoor scenarios. 

. TR beamforming techniques, which have been thoughtfully analyzed in other scenarios 
and take advantage of rich scattering, have not been studied int he context of mm-wave 
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massive MIMO. More specifically, [T35TI and Il29ll suggest that SINR in conventional TR 
grows linearly with the number of antennas, enabling low complexity receivers. 

In this context, the contributions of this paper are the following: 

. We introduce a simple channel model for 60 GHz massive-MIMO, which is based on the 
IEEE 802.llad model EDl . We define the probability distribution of the channel taps, their 
PDP, and spatial correlation. 

. We study the performance of conventional TR in multi-user systems when the number of 
antennas at the transmitter is very large. Moreover, we generalize the ETR |[29l to multi-user 
systems, and compare its performance with conventional TR. We demonstrate that, provided 
a sufficiently large number of transmit antennas, TR does not need further equalization, 
becoming an attractive beamforming alternative. 

. Using the previous analysis, where we find that conventional TR performance is IUI-limited, 
we propose a novel TR multi-user beamforming technique that minimizes IUI and exploits 
rich multipath commonly found in indoor environments. We call this technique interference¬ 
nulling time-reversal (INTR). 

. We analyze and compare numerically the performance of these TR techniques using the 
proposed statistical MIMO channel model for 60 GHz. 

Commonly Used Acronyms in this Paper 

AP - Access Point; CB - cubicle scenario; CIR - channel impulse response; CR - conference 
room scenario; ETR equalized timer-reversal; INTR - interference-nulling time-reversal; ISI - 
inter-symbol interference; IUI - inter-user interference; LR - living room scenario; MPC - mul¬ 
tipath component; PDP - power-delay profile; TR - time-reversal; US - uncorrelated scattering. 

Notation 

Lower and upper case symbols represent signals in the time and frequency domains, re¬ 
spectively. Boldface symbols represent vectors or matrices, whose dimensions are specified 
explicitly, (x) is the convolution operator between two signals. E[ ] represents expectation over a 
random variable. The operators (-) T , (•)*, (-) H and (-) -1 represent transpose, complex conjugate, 
Hermitian transpose, and matrix inverse, respectively. The norm of the vector a is denoted as 
||a|| = y/ (a, a), where (a, b) = b"a represents the complex inner product of vectors a and b. 
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The superscripts (■) tr , {■) eq , and (-) m denote variables calculated using time-reversal, equalized 
time-reversal, and interference-nulling time-reversal pre-filters, respectively. 

II. Time-Reversal Beamforming System Model 

In this section, we present the general discrete signal model for TR beamforming. We first gen¬ 
eralize to the multi-user case two TR techniques for single-user scenarios lf29l . These techniques 
serve as a baseline comparison for the novel INTR introduced in Section [Till 

A. General TR Signal Model 

Consider a digital baseband downlink wireless communication system, consisting of one 
Access Point (AP) with M transmit antennas and N single-antenna user terminals as depicted 
in Fig. Q] The transmitter has a very large number of antennas, so M » N. We denote the 
transmit antenna set as Ad = {1,2,..., M } and the user set as A f = {1,2,..., }. Also, let 

m, m' e A4 and n, n' e Af be arbitrary elements in those sets. The AP transmits simultaneously 
an independent data stream to each user. Let s n (t) be the complex random signal transmitted to 
the n-th user, where t e Z + is the discrete time index. These transmitted signals are assumed 
to have unit average power, i.e. E [|s n (f)| 2 ] = 1, Vn,t, regardless of the modulation. In a TR 
multi-user system, the transmitter sends independent signals simultaneously to the users using 
different pre-filters for each one of them. Thus, the baseband transmitted signal from the m-th 
antenna is 

N 

Xm(t) = VPJ] S n (t) <8>*4 in (-t), (1) 

72=1 

where p is the total average transmitted power in the AP, p m , n (t) is the power-normalized 
pre-filter from the m-th transmit antenna to the n-th user (with a duration of L p samples, i.e. 
t = 0,..., L p — 1), and h m , n {t ) is the random channel impulse response (CIR) from the m-th 
transmit antenna to the n-th user (with a length of L samples). The random CIR vector to the 
n-th user is defined as 

h„(t) = [hi,nit), ..., h M ,n(t)] T e <C M (2) 

In Section [TV] we introduce the statistical characterization of h m , n (t) for mm-wave channels. 
Let H m>n (f ) be the discrete Fourier transform (DFT) of h m , n (t). In an analogous way to the 
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Fig. 1. System model. An AP with M transmit antennas sends simultaneously an independent data stream to N single antenna 
users using time-reversed pre-filters p* „(-(). 


time domain representation, the steering vector to the n-th user is 

H n (/) = [#!,„(/), • ■ • , H M M)f 6 C M . (3) 

The selection of p m ,n(t ) depends on the particular TR technique, as discussed later in this 
section. We define the pre-filter vector to the n-th user as 

P n(t) = [Pl,n(t), ■ ■ ■ ,PM ) n( t )] T 6 C M . (4) 

Let P m ,n(f) t> e the DFT of p miTl (t). Then, we define the frequency domain pre-filter vector to 
the n-th user as 


Pn(/) = [A,n(/),...,PM,n(/)] T eC M . 


( 5 ) 




















The received baseband signal at user n is 


M 


Vn(t) = VP S n{t) ® Y Pm,n(-t) ® h m,n(t) 


771=1 


signal directed to the n-th user 
N M 

d" -\[p ^ 'Yl Sn '(ty ® Pm,n'( — t) ® h"m,n{t) 


n'= 1 m=l 
n'^n 


■v— 

IUI 


+ ^ n (t) 


( 6 ) 


where 2 n (£) represents additive white Gaussian noise (AWGN) with variance cr z . Next, we extend 
the conventional TR and equalized TR single-user formulation in [|29l to multi-user scenarios 
by explicitly defining the pre-filter p m ,n(t ) in terms of the CIR. Note that p m , n {t) is properly 
normalized so the transmitted power p regardless of the number of antennas or users. 


B. Multiuser Conventional TR Beamforming 

The general idea behind TR is to use the time-reversed CIR from every antenna to the receiver 
as a pre-filter for the transmitted signal. Such pre-filter acts as a beamformer in the spatial domain, 
focusing the RF signal around the receiver. In conventional TR, assuming perfect channel state 
information (CSI) at the transmitter, the pre-filter vector is 

h n(L - 1 + t) 


Pn (t) = 




(7) 


where Pjf is a normalization factor introduced to ensure that the total transmitted power remains 
constant in every realization, this is 


N L—l 


p lt - 2 S « h " 


( 8 ) 


n =1 1=0 


Note that, in this case, the pre-filter’s length is equal to the CIR length, i.e. Lf = L. Replacing 
the conventional TR pre-filter into © and using the definitions in Section III-A1 the time domain 








received signal in conventional TR is 


V f-i 


■L-l 


+ 


+ 


1=0 



desired symbol directed to the n-th user 
2L—2 M L-l 


EES h m ,„(l')h' m , n (L -1-1 + l') s „(t - l) 


1=0 m=1 l'=0 
&L -1 


ISI directed to the n -th user 



N M 


Y Y h kn'( L ~ 1 - *) ® h mi n(t) (X) s' n (t) 


n '=1 m= 1 
n'^n 


IUI 


+ -n(^) • 


(9) 


This received signal is composed of four terms: i) the desired symbol multiplied by a real 
factor resulting from coherent combination of multipath components in the CIR, ii) ISI caused 
by incoherent addition of CIR components, Hi) IUI caused by the signals directed to other 
users (whose TR pre-filters do not match the CIR to the n-th user), and iv) AWGN. Thus, in a 
conventional multiuser TR beamforming system, ISI and IUI are important problems that hamper 
detection. In the single-user scenario, ETR was proposed before as a solution to mitigate the ISI 
component in the received signal |[29l . We extend ETR to the multiuser case next. 


C. Multi-user Equalized TR Beamforming 

ETR uses the TR pre-filter in cascade with a ZF pre-equalizer in order to mitigate the ISI of 
conventional TR. In ||29l , it is demonstrated that ETR outperforms conventional TR with respect 
to BER in a single-user scenario, with a marginal loss in the spatial focusing capability. We now 
extend this technique to the multi-user scenario by defining the pre-filter vector components for 
the n-th user as 

h m , n (L — 1 + t) (x) g*[—t) (10) 


PZnit) = 


Vp, 


eq 

h 


where g n (t ) represents a ZF linear equalizer with length L E . Thus, we have Lf = L + L E — 1. 
The normalization factor is 

M N L+Le-2 

E K, n (L-l-l)®gM\ 2 . (11) 


p e 9 _ 
r h — / 1 


m= 1 n=1 1=0 
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One equalizer is required for each user, with the n-th equalizer designed to satisfy 

M 

Quit) <8> 2 h ™,n( L h m , n (t) = 5(t - to), (12) 

m= 1 

where t 0 is an arbitrary delay. Equation (fT2l) can be written as an over-determined system of 
linear equations on g n (t), t = 0,..., L E — 1. Thus, perfect ZF equalization is not possible with 
a finite equalizer’s length 06ll . but a good approximation can be achieved with a sufficiently 
large Le, eliminating the second term in ©. A detailed discussion on this subject is provided in 
|[29l . Using the ETR pilot, and assuming perfect equalization, the time domain received signal 
at user n is 

y*{t) = s *(* ~ *o) 

'---V-' 

signal directed to the n-th user 

I - N M 

+ \~p^ 2 2 s n'(t)®gn'(t) (X) h* m y(L - 1 - t) (X) h miU (t) 

V h n / =1 m= i 

n' j=n 

" -v-" 

IUI 

+ z n (t) , (13) 

noise 

which has no ISI term, but still contains IUI. Thus, both conventional TR and ETR performance 
is limited by ISI and/or IUI, as detailed next. 


D. Performance Analysis of TR and ETR 

We now turn our attention to the power components in ([9]) and (fl3l) . following the same pro¬ 
cedure as in [|29fl . The fundamental assumptions are that the system operates under uncorrelated 
scattering with uncorrelated channels between users, and normalized channel power, as stated in 
Section UV] In the derivations below, we also employ the approximation E[«/6| E[a]/E[6] for 

two random variables a and b, as analyzed in Il29l . 051 for TR systems. Complete derivations are 
not shown due to space constraints. Here, we verify that TR is a suitable technique for massive 
MIMO systems, where M » 1. Let Pf r , P-f, and P-f represent the power in the first, second, 
and third terms in ©, respectively. Then, the average desired signal power is 

L—l 


E [Pf] = E 


P 

r>tr 

l h. 


2 ii h - 


s n (t — L + 1) 


1=0 


MpT 

N 


( 14 ) 
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The average ISI power in © can be approximated as 


E [f£] = E 


TDtr 

r h 


2L-2 M L—l 


2 2 2 h m Aiy cal - 1-1 + r) Sn (t -1) 


1=0 m= 1 l'=0 
1+tL -1 


P 


2L-2 M M L—l 


MNT 


2 2 2 2 E [ ft -».»( ( >m>(0i] 


Z=0 m= 1 m'= 1 Z'=0 
l^L-1 


E [h* m>n (L -1-1 + l')h m ,, n (L -1-1 + /')] . 


(15) 


An approximation to the average IUI power is 




N M 

2- 

e rej= e 

p 

Pp 

2 2 h *m,n'{ L - 1 - *) ® ® s' n (t) 



1 h 

n '=1 m =1 
n'^n 



P 


N 2L-2 M M L—l 


MNT 


2 2 2 2 2 E [a™, ++L ,„(!')] 


n'= 1 1=0 m= 1 m'= 1 l'=0 
n'^n 


xE [h* m>n ,{L -1-1 + l')h m ,, n ,(L -1-1 + 1')] . 


(16) 


For ETR, the average desired signal power is bounded by 

E [P* q ] = E 

The average IUI power in ETR has a similar form to (fT6l) . but it is not shown here since it is not 
the focus of this work. However, we analyze it numerically in Section [Vj From (fl4l)- (fl6l) . we 
can make the following remarks with respect to TR beamforming in massive MIMO systems: 

* Both ISI and IUI powers are highly dependent on the propagation conditions. More specifi¬ 
cally, power delay profiles and correlation between antennas are present in the terms of the 
form E [hm^nipjh*^ n {l)]. Thus, increasing spatial correlation would increase both ISI and 
IUI, degrading performance. 

* In the case of uncorrelated antennas, E n {l)] = 0 if rn + m!. Hence, the sums 
would only depend on the power delay profile (which is the same for all antennas), and 
both ISI and IUI powers would be independent of M. 

* Desired signal power increases linearly with M. Thus, in uncorrelated channels E- 
oo and E[Pf /P^] —> oo as M —> oo. This implies that, with a sufficiently large number of 


< 


MpT 

N 


(17) 
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antennas, a conventional TR beamforming system is noise limited instead of interference 
limited. 

* However, if channels are spatially correlated (as in realistic scenarios), equalization and 
interference mitigation provided by TR reduce. 

* Note that, given that CIR statistics for users n and n' are the same, i.e. E \h* n n (l)h m > >n (l)] = 
E \h* m n '(l)h rn i. n '(l)] V/, then PP is larger than PP by a factor on the order of the number 
of users. Thus, IUI mitigation should be given priority over equalization when proposing 
improvements over conventional TR. 

Given this characteristics of TR beamforming in massive MIMO, we now propose a novel 
TR extension to overcome the problems of IUI, even under highly correlated channels. 

III. Interference-Nulling Time-Reversal Beamforming 

We are now concerned with the design of pre-filter vectors that combine the spatial focusing 
properties of conventional TR, while also providing additional IUI mitigation. We start from the 
frequency representation of the received signal, and formulate an optimization problem for the 
design of the pre-filters. The frequency domain equivalent of © is 

Y n (f) = VP <H fl (/),P n (/)>S n (/) 

"■-V-' 

signal directed to the n-th user 
N 

+ Vp2< h »(/).PA/)>SW(/) 

n '=1 
n! j=n 

" -V-' 

IUI 

+ Zntn, (18) 

noise 

where S n (f) is the DFT of s n (t), and Z n (f ) is the DFT of z n (t). Appropriate zero padding is used 
in the time domain in order to represent linear convolution as a product in the frequency domain. 
The complex inner product defined above allows a convenient simplification in (fl8l) with respect 
to (©, which is useful for the problem formulation. Let H(/) = [Hi(/)... Hjv(/)] e C MxJV be 
the matrix with columns given by the steering vectors to all users. Also, let H_„(/) e C Mx/Y 1 
be the matrix formed by removing the n-th column from H(/), i.e. removing the steering vector 
to user n. For notational simplicity, we drop the frequency dependence in the remainder of 
this section. Note that the IUI power in (fl8l) is proportional to Yu n '*n | (H n , P n /) | 2 . Thus, our 
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Fig. 2. Geometrical interpretation of the optimization procedure. The optimum pilot in the frequency domain is the conventional 
TR prefilter projection onto the nullspace of H? n . This ensures that IUI is set to zero for every user and every frequency. 


objective is to find the pre-filter P* n which is closest to the conventional TR solution in the 
frequency domain (providing partial equalization of the received signal), and such that the IUI 
is set to zero. Formally, this optimization problem can be formulated as 

P*,n = argmin ||P£ - P n f 

t n 

subject to H ff n P n = 0, 

\/neAf , V/ e [0,..., L + L p - 1], (19) 

whose solution is 

P*, n = (l M ^ H_ n (H^FU)' 1 H" ) P* 

Vn,/. (20) 

where I m is the M x M identity matrix. Thus, we call P* „ the interference-nulling time-reversal 
(INTR) pre-filter in the frequency domain. Geometrically, the constraint in the problem ensures 
that the vector P* n e null { H u _ n } , and the solution is the projection of P^' into that null space. 
This is illustrated in Fig. (2j 

IV. Channel Model for 60 GHz Massive MIMO 

As mentioned above, TR actually benefits from rich scattering, so it can be conveniently 
applied for indoor wireless communications. In this section, we briefly describe the IEEE 
802.1 lad model for 60 GHz SISO systems in such scenarios lf2TI . and extend it to the correlated 
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multi-user massive MIMO case. In the following, we use a statistical description of h m<n (t), 
given by its probability distribution, power delay profile (PDP), and spatial correlation in the 
context of massive MIMO systems. 

A. Channel Tap Distribution 

The most popular channel model for mm-wave propagation is the IEEE 802. llad. This is 
a SISO double directional statistical channel model based on a limited set of measurements, 
complemented with ray-tracing simulations. This model is defined for three indoor scenarios: 
conference room (CR), living room (LR), and cubicle environment (CB). Some important model 
features include: support of two types of antennas (isotropic and basic steerable antenna array), 
support of polarization, wideband and pathloss modeling under LoS and NLoS situations. 

The IEEE 802.1 lad channel model follows a scattering cluster structure, both in time and 
angular domains. Thus, several multipath components (MPCs) observed in the CIR have similar 
propagation delays and angles of departure/arrival. More specifically, each central ray arriving at 
the receiver has pre-cursor rays (which arrive earlier) and post-cursor rays (which arrive later). 
This is due to irregular scattering objects and geometrical features which are large compared to 
the wavelength. Both pre-cursor and post-cursor rays have less amplitude than the central ray. 
The resolvability of those MPCs depend exclusively on the system’s sampling time (bandwidth). 
When those MPCs are not resolvable, they contribute to the same tap in the CIR. Given this 
propagation characteristics, we assume that h m , n (t) has zero mean and that \h m ,n(t)\ is Nakagami 
distributed, with parameters m and Q, for all m, n and t [f37Tl . Recall that the m parameter in 
the Nakagami distribution is analogous to the K factor in the Rician distribution, and that a 
larger m implies a large power ratio between the central ray (specular component) and the other 
rays (diffuse components). The parameter Q depends also on the amplitudes of the specular and 
diffuse components, and on the channel PDP (tap average power) If37l . Table U shows the values 
of m and RMS delay spread in the IEEE 802.llad scenarios. The larger value of m in the CB 
scenario is due to the reduced scattering within the cubicles, which reduces the number and 
power of diffuse components contributing to each channel tap. 
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TABLE I 

NAKAGAMI m PARAMETER AND RMS DELAY SPREAD OF IEEE 802.11 AD SCENARIOS 


Scenario 

Nakagami m parameter 

RMS delay spread [ns] 

CB 

4.34 

3.47 

CR 

2.56 

4.82 

LR 

1.74 

7.81 


B. Power Delay Profile 

We are particularly interested in the PDP, a second order statistic defined as 

A h (t) = E[\h m}n (t)\ 2 ], Vm, n, (21) 

where the expectation is calculated over CIRs that are subject to the same large-scale fading 
lf38l . Signal power components depend on the PDP and spatial correlation, as seen in Section 
III-DI We assume that all CIR in the system have the same PDP This is valid for mm-wave 
indoor environments, where APs are usually positioned on or close to the ceiling and similar 
shadowing affects all elements in the transmit array. We also define the following constraint on 
the CIR total power: 

L—l L—l 

^EjJVnWI 2 ] = 2^^ =r ’ (22) 

4=0 4=0 

where T « 1 is a constant accounting for channel induced propagation losses. This constraint 
implies that all channels between the transmit antennas and each receiver have same average 
power. Fig. [3] shows PDPs obtained over 10 6 realizations of the IEEE 802.llad model for the 
three scenarios simulated under NLoS and isotropic antennas. We do not consider LoS situations 
since they correspond to flat-fading channels, which are of no interest here. Isotropic antennas are 
assumed so the system can take advantage of all MPCs in the channel. In practice, planar omni¬ 
directional antennas (e.g. 11391 ) would be a good alternative for implementation. We observe that 
RMS delay spread is minimum for the CB scenario, where an AP is located in the ceiling of an 
office populated with cubicles. In that case, scattering is confined within the cubicle’s structure 
and other delayed paths (e.g. reflections from outer walls) are obstructed. On the other hand, 
CR and LR scenarios correspond to more open spaces, where first and second order reflections 
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from walls are considered. Those reflections cause long tails in their PDPs, increasing their delay 
spread. 

C. Spatial Correlation Model 

Consider m, m! e Ai, n, n' e J\f, and t, t' e {0,. .., L— 1}. We make the following assumptions 
with respect to CIRs in the systems: 

. CIR are correlated across transmit antennas, i.e. the spatial channel autocorrelation function 
is Rh(Ad) A 0, where Ad is the distance between two measured CIRs. The specific 
correlation structure depends on the array configuration, but it is assumed that the process 
in wide sense stationary with respect to the space. This implies that E n (£)] A 0. 

. Different users have uncorrelated CIRs to the AP, i.e. h myn (t) and h rn n i(t) are uncorrelated 
if n A n!, Vm, t. This is due to the fact that MPCs are independent for different users. This 
can be clearly seen in the CB environment, where each user is assumed to be in its own 
cubicle. 

. CIR taps are uncorrelated, i.e. h m ,n{t ) and h m , n {t') are uncorrelated if t A t 1 , \/m,n. This 
is the conventional uncorrelated scattering (US) assumption widely used in the literature 
ED, and implies that contributions to different taps come from different scatterers. 

Nakagami correlated variables (across antennas) are generated according to the method de¬ 
scribed in lf37l . as follows. Consider the setting in Fig. BJ A planar randomly-oriented array 
with M isotropic elements is located in the environment according to the standard lITill . An 
isotropic receiving antenna is randomly located in the environment as well. Each tap is assumed 
to have specular and diffuse contributions from an irregular scatterer (located according to 
the corresponding delay), whose amplitudes depend on the PDP and the desired m parameter. 
All contributions to a fixed tap in a given CIR come from the same scatterer, with different 
taps corresponding to different scatterers. Using this procedure, the resulting normalized spatial 
correlation function R h (Ad) is shown in Fig. [5} These results are consistent with measured and 
simulated spatial correlations in 60 GHz channels, e.g. |[22l . High correlation values are caused 
by the reduced number of dominant MPC contributing to each tap. The specific correlation 
between transmit antenna elements depends only on the geometry of the array. For the numerical 
validation shown in Section [V] we use rectangular arrays with 32 (8 x 4), 64 (8 x 8), or 128 
(16 x 8) elements with a uniform separation of 20 mm. 
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Fig. 3. Power delay profile of IEEE 802.llad channel model scenarios with isotropic antennas. RMS delay spreads are 3.47 
ns for the CB scenario, 4.82 ns for CR and 7.81 ns for LR. 



Fig. 4. Method to generate correlated Nakagami CIR. Different taps are assumed to have contributions from specular and 
diffuse reflections from different objects. The transmit array is planar (rectangular) with uniformly distributed elements. 


V. Numerical Results and Discussion 

In this section, we present numerical results for the performance analysis of the multiuser TR, 
ETR and INTR techniques, as described in Sections [TT] and Hill 

A. Pilot Length and Channel Correlation 

First, we analyze the impact of pre-filter’s length L p and spatial correlation on the signal 
power components. We calculate the values of P s , Pi S i, and P lul for the three techniques over 
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Fig. 5. CIR spatial correlation as a function of distance Ad, calculated over 1000 CIR realizations. 

1000 channel realizations with and without spatial correlation. Results are shown in Table [IT] for 
2 and 10 users in the CB scenario with 64 antennas. For conventional TR beamforming, it is 
clear that IUI power is the main problem for multi-user communications, as it can be up to an 
order of magnitude greater than ISI power. It is also observed that channel correlation decreases 
desired signal power and increases interference, affecting the overall system performance. Also, 
note that when using spatially uncorrelated channels, both ISI and IUI suffer small or no change 
when increasing the number of antennas. However, when the CIRs are correlated both types 
of interference suffer a small increase. For ETR, increasing pre-filter’s length improves ISI 
suppression, but IUI remains the same as in conventional TR. ISI reduction in ETR by increasing 
L p is a typical consequence of zero-forcing equalization ll29l . However, ETR is not designed to 
mitigate IUI. Thus, BER performance of TR and ETR are expected to be very similar since the 
scenarios we consider are clearly IUI limited. 

For INTR, IUI mitigation improves by increasing L p . This is due to the discarding of L — 1 
time samples when performing the transformation between the frequency domain prefilter (of 
length L + L p — 1) and the time domain prefilter (of length L p ). Such discarding is necessary due 
to the circular convolution theorem. Thus, the time domain prefilter is a least squares projection 
of the optimum frequency domain solution. The error in the projection reduces as L p increases. 

We observe the impact of signal power components over the BER performance in Fig. [6] 
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Fig. 6. Prefilter length ( L p ) vs (a) ISI power in ETR, and (b) IUI power in INTR. These results were obtained with L = 60 
and N = 10 in the CB scenario. Other signal components in each technique remained approximately constant vs. L p . It is 
noted that increasing L p reduces ISI power in ETR and IUI in INTR. This is due to the discarding of L — 1 time samples when 
performing the transformation between the frequency domain prefilter (of length L + L p — 1) and the time domain prefilter 
(of length L p ). Such discarding is necessary due to the circular convolution theorem. Thus, the time domain prefilter is a least 
squares projection of the optimum frequency domain solution. The error in the projection reduces as L p increases. 


where the influence of channel spatial correlation is also shown. Signal to noise ratio is defined 
as SNR = fN/o'j, where ai is the variance of z n (t ) Vn, t. These results were obtained for 5 users 
and 32 antennas in the CB scenario, with a transmission of 10 6 BPSK symbols over 1000 channel 
realizations. Performance of both TR and ETR is limited by IUI, which causes a lower bound 
in the BER. We notice that ETR does not provide a significant improvement over conventional 
TR in the case of multi-user massive MIMO systems. Thus, ETR does not offer any advantage 
for such scenarios, given its greater computational complexity with respect to TR. In the case 
of INTR, IUI is successfully mitigated and hence INTR outperforms the other techniques. We 
also observe that channel correlation degrades system performance in all cases. 

B. Number of Antennas and Number of Users 

Fig. [8] shows the average BER performance results per user for varying number of antennas 
and users. These results where obtained with the transmission of 10 6 BPSK symbols over 1000 
spatially correlated channel realizations. We used a fixed pre-filter length L p = 90 and the 
CB scenario PDP Conventional TR performance results are consistent with the analysis made 












20 


TABLE II 

Received signal power components for M = 64. Values are normalized /(pr). 


Number of users N 

Technique 

Lp 

— L = 

60 

L p = 90 

L p = 120 

Ps 

Pisi 

Piui 

Ps 

Pisi 

Piui 

Ps 

Pisi 

Piui 

Uncorrelated channel 


TR 

32 

0.15 

0.51 

- 

- 

- 

- 

- 

- 

2 

ETR 

- 

- 

- 

31.9 

0.01 

0.52 

31.9 

0.001 

0.52 


INTR 

31.7 

0.15 

0.09 

31.6 

0.16 

0.01 

31.6 

0.15 

0.002 


TR 

6.4 

0.03 

0.9 

- 

- 

- 

- 

- 

- 

10 

ETR 

- 

- 

- 

6.37 

0.002 

0.9 

6.37 

0.0003 

0.9 


INTR 

5.77 

0.04 

0.17 

5.58 

0.04 

0.02 

5.55 

0.04 

0.004 

Correlated channel 


TR 

31.4 

0.6 

0.88 

- 

- 

- 

- 

- 

- 

2 

ETR 

- 

- 

- 

30.3 

0.02 

0.88 

30.3 

0.003 

0.88 


INTR 

30.8 

0.61 

0.15 

30.6 

0.6 

0.02 

30.6 

0.6 

0.003 


TR 

6.39 

0.15 

1.48 

- 

- 

- 

- 

- 

- 

10 

ETR 

- 

- 

- 

6.07 

0.004 

1.47 

6.07 

0.0005 

1.47 


INTR 

5.32 

0.14 

0.25 

5.06 

0.14 

0.04 

5.02 

0.14 

0.006 


in Section III-Dl The desired signal power increases linearly with the number of antennas while 
interference components remain constant. Thus, the minimum achievable BER per user improves 
by increasing the number of antennas, providing diversity gain. On the other hand, increasing 
the number of users with a fixed number of antennas decreases the desired signal power and 
increases IUI. This is reflected in a higher BER for larger N. INTR outperforms conventional 
TR in every simulated scenario. Nevertheless, the performance improvement provided by INTR 
is more evident with a large number of users or a limited number of antennas. 


C. Average Achievable Sum Rate 

The achievable sum rate measures the downlink spectral efficiency in a multiple-access system. 
Assuming that each user treats ISI and IUI as Gaussian interferences, and according to the model 
defined in Sections [TT] and [Oil the average achievable sum rate for a multi-user TR system is 


N 


2 l0 §2 1 + 


n=l 


s,n 


p. . _l p. . _L 

1 isi,n ' 1 im,n ' u z 


R = ]E 


(23) 
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w 



SNR [dB] 


Fig. 7. Average BER per user comparison of TR, ETR, and INTR, under correlated and uncorrelated channels (across antenna 
elements) with M = 32, and N = 5. Results are shown for different pilot lengths ( L p ). It is observed that spatial correlation 
increases ISI and IUI, degrading performance. Also, increasing prefilter’s length improves IUI mitigation in INTR. 


where P StU , P lsl . n , and P, w , n are the desired signal power, ISI power, and IUI power, respectively, 
calculated at user n for a given realization. For simplicity, it is also assumed that the channel is 
used for downlink transmission all the time. A proper reduction factor can be used to account 
for uplink time in a TDD system or channel estimation overheads. Numerical results for the 
average achievable sum rate are shown in Fig. [9l These results were obtained in the LR scenario 
with correlated channels and L p = 90. As seen, INTR offers a significant improvement over 
conventional TR, doubling its rate in some cases and providing a remarkable multiplexing gain. 
In addition, we simulated a more extreme case with N = 30, and N = 50 and 128 antennas, 
with the purpose of further demonstrate the capabilities of TR to handle IUI. Results are shown 
Fig. [9c] Even though the assumption of uncorrelated CIR between users is hardly met when 
N is that large, results show that an outstanding efficiency of more than 170 bps/Hz can be 
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(a) (b) 


Fig. 8. Avera BER per user for TR and INTR. (a) Different number of antennas M with L p — 90 and N = 5. (b) Different 
number of users N with L v = 90 and M = 64. An important diversity gain is achieved even in spatially correlated channels. 
The effect of IUI is mitigated by increasing the number of antennas. 
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SNR [dB] 
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(C) 


Fig. 9. Achievable rate of TR and INTR in the LR scenario, (a) M = 32 antennas, (b) M = 128 antennas, (c) M = 128 
antennas with an extreme number of users. The multiplexing gain increases with the number of antennas. 


achieved with INTR. In all the simulated scenarios our proposed INTR technique outperforms 
conventional TR, as it can better withstand an increase in user load. 

VI. Conclusion 

We have analyzed a baseband TR beamforming system for mm-wave multi-user massive 
MIMO. We studied conventional TR and equalized TR and found that their performance is 
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IUI limited. We also noticed that, when the number of antennas is large, the ratio between the 
desired signal power and ISI or IUI power increases. Thus, we confirm the potential of TR as 
a beamforming technology for massive MIMO. We also note that equalizing solutions such as 
ETR are not necessary when the number of transmit antennas is large. After identifying IUI as 
the main detection impairment for TR systems, we propose a modified technique called INTR. 
This technique calculates the transmit pre-filters in the frequency domain that set the IUI to zero 
and are closest to the original TR solution. We proposed a 60 GHz MIMO channel model, where 
CIR taps are modeled with Nakagami distributed amplitudes. In addition, we use PDPs given by 
the IEEE 802.1 lad SISO NLoS model, and generate spatial correlation in the CIRs according to 
a geometrical model. By means of numerical simulations, we verified that the proposed INTR 
outperforms conventional TR with respect to average BER and achievable sum rate. In particular, 
we note that INTR performance is extremely tolerant to increases in the number of users, and 
provides both diversity and multiplexing gains simultaneously. 
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